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THEORY  OF  MULTICAVITY  GYROKLYSTRON  AMPLIFIER 
BASED  ON  A  GREEN’S  FUNCTION  APPROACH 


INTRODUCTION 


A 


_ „  jyroklystron  amplifiers  show  great  promise  as  high-power  and  high-gain  devices  in  the 
microwave  and  millimeter  wavelengths,  ''lory  et  al.  [1]  performed  the  first  gyroklystron  amplifier 
experiment  in  a  two-cavity  configuration  with  cylindrical  TE°m/  TE°m  modes.  Recently,  Bollen  et 
al.  [2]  have  performed  a  three-cavity  gyroklystron  amplifier  experiment  operating  with  rectangular 
TE’Joi  mode  at  4.2  GHz.  They  achieved  an  output  power  of  54  kW  with  16%  efficiency  and  0.2% 
bandwidth.  To  achieve  higher  gain  and  moderate  bandwidth,  a  multicavity  configuration  with  stag¬ 
gered  tuning  will  be  necessary  as  in  conventional  klystrons.  The  design  of  such  a  configuration 
requires  the  optimization  of  many  parameters.  The  large  signal  analysis  and  numerical  codes  [3-5] 
developed  previously  for  two-cavity  gyroklystrons  can  be  extended  to  the  multicavity  configuration, 
but  the  procedure  would  require  extremely  large  computation  time.  Moreover,  the  bandwidth  and 
the  effect  of  staggered  tuning  could  not  be  calculated  since  the  radio  frequency  (RF)  fields  are  not 
determined  self-consistently.  Therefore,  it  is  necessary  to  develop  an  analytical  linear  theory  for 
multicavity  gyroklystrons.  The  linear  theory  will  be  applicable  to  the  prebunching  cavities  where 
RF  field  is  small,  and  these  results  may  be  used  as  an  input  for  large  signal  calculation  in  the  power 
extraction  cavity  to  obtain  saturation  gain  and  efficiency.  Symons  and  Jory  [6]  developed  a  small 
signal  theory  for  the  two-cavity  gyroklystron  by  using  a  sinusoidal  RF  profile  and  obtained 
equivalent  circuit  parameters.  Caplan  [7]  outlined  a  small  signal  self-consistent  theory  of  gyrokly¬ 
stron  amplifier  based  on  Maxwell- Vlasov  theory.  In  this  report,  we  develop  a  self-consistent  theory 
of  the  multicavity  gyroklystron  amplifier  in  which  the  electron  motion  is  represented  by  the  gen¬ 
eralized  pendulum  equation  [8,9].  The  derivation  of  these  equations  [10]  are  outlined  in  the  next 
section.  These  equations  are  solved  in  the  third  section  for  the  multicavity  gyroklystron  configura¬ 
tion  by  using  a  Green’s  function  approach  to  satisfy  appropriate  RF  field  boundary  conditions. 
Small  signal  results  are  then  obtained  by  the  method  of  successive  approximation.  Effects  of  axial 
beam  velocity  spread  are  also  included.  In  the  fourth  section  we  calculate  the  small  signal  perfor¬ 
mance  characteristics  of  the  three-cavity  gyroklystron  amplifier  configuration  used  by  Bollen  et  al. 
[2], 


SELF-CONSISTENT  EQUATIONS  FOR  GYROTRONS 

Figure  1  illustrates  the  multicavity  gyroklystron  amplifier  under  study.  The  cavities,  separated 
by  drift  regions,  are  placed  in  a  uniform  magnetic  field  B0  applied  parallel  to  the  axes  of  the  cavities 
(z-axis).  Let  rwj  and  Lj  denote  the  radius  and  the  length  of  the  jth  cavity;  dj  is  the  length  of  the 
yth  drift  region.  The  entrance  of  the  jth  cavity  is  at 

z  -  z;  -  +  d,) . 

An  input  RF  signal  and  a  monoenergetic  beam  of  electrons  are  injected  from  the  left  into  the  first 
cavity.  The  electrons  follow  helical  trajectories  due  to  the  strong  uniform  magnetic  field  and 
experience  perturbing  RF  fields  in  the  cavities.  It  is  assumed  that  the  RF  fields  are  completely  cut 
off  in  the  drift  regions.  Space  charge  effects  are  also  neglected. 
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Fit-  1  —  A  three-cavity  gyroklystron 


Consider  the  beam  interaction  with  the  electric  field  of  a  TE-type  circular  waveguide  mode. 
Under  most  gyroklystron  operating  conditions  (electron  velocity  vt  much  less  than  phase  velocity 
\p) ,  the  electron  beam  interaction  with  RF  magnetic  field  may  be  neglected.  For  propagation  of  a 
single  TEm,  mode,  the  electric  field  is  given  by 


E,  -  RzlCmn[kHmJ'm(k,mr)e6  +  (im/r)/m(*Mr)e,}F(z)e'<-'-"*>  1,  (1) 

where  <o  is  the  wave  frequency,  Jm  is  a  Bessel  function  of  order  m,  the  prime  denotes  differentia¬ 
tion,  and  km,,  is  the  transverse  wave  number.  The  normalization  constant  is  given  by 

Cnw  “  lUGd,  -  m!)ll,!4(jtm)]'1,  (2) 


Xmn/r*, 


where  xM  is  the  nth  zero  of  J'm.  The  axial  dependence  of  the  RF  fields  in  the  cavities  is  given  by 
the  complex  profile  function 


F(z)  -  |F(z)|*-'f(*>.  (4) 

d\F\ 

|F(z)|  and  dij(z)/dz  are  assumed  to  be  slowly  varying  functions  of  z  such  that  Xf— — ■  <  |F|  and 

dz 

similarly  for  the  phase.  Here  A,.  —  2 ir  v2/«o. 


In  uniform  external  magnetic  field,  the  axial  momentum  ( pt )  of  the  electrons  remains  con¬ 
stant  if  the  RF  magnetic  field  is  neglected.  For  small  values  of  y  —  [1  +  pi  +  pz2]l/2,  we  will 
assume  that  v2  is  approximately  constant.  For  operation  near  cyclotron  resonance  (i.e., 

jQ0 

<o - «  w),  it  is  convenient  to  write  the  transverse  components  ipx,  py)  of  the  electrons  in 


the  following  form: 


px  -  -Pi(t)  sin  +  0 (r)}, 

Py  “  Pi  0)  cos  {fl  U)t  +  0(f)}, 


where 


fl  —  \e\Bofntoy  -  do/y. 


t  ~  to. 


and  fl ,  Pi  ,  and  <t>  are  slow-time  scale  variables;  t0  is  the  time  the  electrons  enter  the  interaction 
region.  Under  these  approximations,  it  has  been  shown  [10,11]  that  the  Lorentz  force  equations 
for  the  electrons  can  be  expressed  approximately  in  the  following  form  in  terms  of  two  slow-time 
scale  variables  pL  and  the  phase  angle  A-wt  -  dr  -  <f> : 


-  -\e\CmlkmlJm-,(kmlro)  J',(kmlrL)\F\  cos  'I', 


I 


m 


t  a  ■ 
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\e\Cmlkm,Jm-s(kmir0)  lf  I  sin  *  +  <»  ~  sfl)pL, 


(10) 


where  s  is  the  cyclotron  harmonic  number  and  r0  and  rL  are  the  guiding  center  radius  and  the  Lar- 
mor  radius.  The  phase  ¥  in  Eqs.  (9)  and  (10)  is  given  by 


V  = 


Sfif 


01  — 


|(f  -  t0)  +  u>t0  -  (  -  <t>  -  (m  -  s)d0 

=  A  —  g  —  (m  —  s)9  o, 
where  0O  is  the  polar  angle  of  the  guiding  center. 

The  wave  equation  for  the  electric  field  profile  function  F(z )  may  be  written  as  [10] 


Crnn  k mn  Jm-Mmnra)<PiJs-\{kmnrL)e  *A>, 


(ID 


4  +  4 

1  - 

/ 

—  fc* 

dz2  c 2 

Q 

n mn 

(16JI0 


.  „  .  (12) 

'woyvMv 

where  I0  =  nol^lvn,av  is  the  beam  current,  and  Q  represents  losses  from  the  cavity.  In  Eq.  (12), 
<0>  denotes  an  average  over  initial  electron  phase  A0  =  oit0,  the  initial  guiding  center  distribution 
and  initial  electron  velocity  distribution  function  /(v): 

<0>  -  f/ms?  f*  X”  OM0M0. 

For  a  "cold"  beam,  <0>  involves  an  average  only  over  initial  phase  and  guiding  center  position. 
For  this  case  vM  av  =  vN. 

For  fundamental  cyclotron  harmonic  operation  (s  =  1),  Eqs.  (9)  to  (12)  can  be  further  sim¬ 
plified  by  noting  that  kmnrL  «  1.  Hence  J[(kmnrL )  and  J\(kmnrL)/ kmnrL  can  be  replaced  by  the 
leading  term  of  the  small  argument  expansion  of  the  Bessel  function;  i.e., 

Ji(kmn //j  =  J\  (kmnti)/  kmnr^  1/2. 

It  is  now  convenient  to  introduce  a  complex  slow-time  scale  transverse  momentum  given  by 

-/A 


P  =  Pie 


03) 


Equations  (9),  (10),  and  (12)  become 


+  _L 
dz  v„ 


Of 


oi  — 


,  /,  \  r-  ~ dm— D®n 

= - =- - Jm-\(kmnr0)  Fe 

^vii 


(14) 


and 


4 

1- 

a>2 

i 

_  l2 

I  dz2  T 

c2 

I'l 

*mn 

Cmn kmn  Jm-\( kmn rQ) < pe  °>  . 


(15) 


y  m0\IKav 
:d  the  rela 

increased  by  a  normalization  scheme  in  which  the  wavelength  Xq  of  the  radiation  is  an  arbitrary 


In  Eq.  (14),  we  have  used  the  relation  -7-  =  vz~.  The  generality  of  Eqs.  (14)  and  (15)  can  be 

dt  dz 
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number.  For  finite  bandwidth  operation,  Ao  will  refer  to  the  center  frequency  of  the  band.  The 
normalization  scheme  can  be  achieved  as  follows  (the  normalized  quantities  are  denoted  by  a  bar): 

(a)  length  normalized  to  A0(z  -  z/A0) 

(b)  velocity  normalized  to  c(/3  -  v/c) 

(c)  frequency  normalized  to  —  <3  - 

M  c  } 

(d)  electric  and  magnetic  fields  normalized  to 

(E  =  eE\o/m0c2,  B0  -  e8oA</ m0c). 

Other  quantities  such  as  kmn,  p,  t  ,  and  F  are  normalized  consistently  to  km„  -  kmn  Ao,  p  -  p/m0c , 
t  —  tcf  A0  ,  and  F  =  eF/m0c2.  After  these  procedures,  Eqs.  (14)  and  (15)  become 


dp  .  i  _  f)0  -  ^n.ov  = 

^  —  p - —E, 


dz  An 

f  d2 


+  k2\E  -  —Uq<p>, 


where 


A  -  —  1  -  -i-  -  A2  An 
2  1  a/wi 


E  2  ^ 


7  _1_  .  g/0  .  .  ^-1  (19) 

0  2ir  e0m0YC3  A?i<JV  (x^„  -  m2)/^(xmn) 

Equations  (16)  and  (17)  constitute  a  set  of  nonlinear  coupled  equations  for  the  gyrotron.  In  the 
next  section  we  give  an  analytical  solution  of  these  equations  in  the  linear  approximation. 

GREEN’S  FUNCTION  APPROACH  AND  LINEAR  THEORY 

We  solve  Eqs.  (16)  and  (17)  for  a  "cold"  beam  where  Ai.ov  “  An-  At  the  end  of  this  section 
we  point  out  the  modification  necessary  to  consider  axial  velocity  spread  of  the  beam. 

The  general  solution  of  Eqs.  (16)  and  (17)  may  be  written  as 


where 


rz—  if-  b(z)dz" 

p(z)  **  Ip(zj)  -  f  F{z')e  zj  d z']e  *) 

*j 

E(z)  —  -il0  fQ  J  G{z,  z')  <p(z')>dz' 
+  Aj  sin  kzj(z  —  Zj  —  Lj), 


Zj  “  £  (L/  +  </(), 
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and  the  detuning  parameter 


A(z) 


JM 

y(z)y 


(22) 


In  Eq.  (21)  the  first  term  is  the  particular  integral  giving  the  contribution  to  the  field  due  to  the 
perturbed  current  density,  and  the  second  term  is  the  homogeneous  solution.  G(z ,  z')  is  the 
Green’s  function  and  should  be  chosen  to  satisfy  the  correct  boundary  conditions.  If  we  assume 
that  E(z)  vanishes  at  both  z  —  zj  and  z  —  lj  +  Lj,  then  G(z,  z')  is  represented  by 


G(z,  z')  -  ^ 


kzj  sin  kzjLj 


(23) 


This  form  of  G(z,  z')  applies  to  all  the  cavities  except  the  first  and  the  last  ones.  In  those  two  cav¬ 
ities,  boundary  conditions  corresponding  to  the  incoming  and  outgoing  waves  might  be  included. 
In  this  report  we  assume  that  Eq.  (23)  applies  to  the  last  cavity  also.  In  the  first  cavity  the  input 
signal  produces  energy  modulation  on  the  electron  beam.  If  the  beam  current  is  not  too  high  and 
the  cavity  is  short,  then  the  change  in  energy  of  the  electrons  in  the  first  cavity  is  small  and  the 
input  signal  is  perturbed  only  slightly  by  the  electron  beam.  Under  these  circumstances,  the  electric 
field  in  the  cavity  is  well  represented  by  the  homogeneous  solution 

E(z)  -  E0 1  sin  kz<i(Li  -  z),  (24) 

where  E0 ,  is  related  to  the  input  power  incident  at  z  -  0.  In  all  other  cavities,  the  RF  field  is 
induced  by  the  bunched  beam,  and  the  homogeneous_part  of  the  solution  in  Eq.  (21)  will  be 
•neglected.  We  have  already  mentioned  that  E(z)  =  0  in  the  drift  regions,  i.e., 
Zj  +  Lj  <  z  <  Zj  +  Lj  +  dj. 


In  the  small  signal  regime,  the  two  integral  Eqs.  (20)  and  (21)  may  be  solved  by  the  method 
of  successive  approximation  where  £01  will  be  considered  a  small  quantity.  The  initial  condition  are 
chosen  to  be 

p( 0)  -  Plo* *  o  <  4> o  <  27t, 

1(0)  -  £oi  sin  (25) 


By  retaining  terms  up  to  first  order  in  E0\,  we  find  from  Eqs.  (20)  and  (24)  that  the 
transverse  momentum  at  the  end  of  the  first  cavity  (z  -  Lx)  of  an  electron  with  initial  phase  <J>0  is 
given  by 


Jl  -  Is5- 

Plo 


(26) 


where 


A0 

X(i) 


JqL]  sin  kzA(Lx  -  D  e^dz' 

£?j{cos  (ki^Lj)  -  cos  (ApI,)}  -1-  )'(A0  sin  kt^L\  - 
— — 


kz  l  sin  A0L|) 


Dp 

v(0) 


J*0  A(z')df' =  Apr  + 


Up 

TO0II 


(27) 

(28) 
(29) 
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8y  _  \p(z)\2-p?0  _  \p(Zj)\2-p?0  p*(zj)  ,Anr 


2y02 


Ki 


i32  1  - 


i 

Qj 


mn 
r2  . 

rWJ 


2yo 

1/2 


VO 


L2E(z')t 

Jzj 


dz\ 


(30) 


(31) 


From  Eqs.  (24),  (29),  and  (30),  we  have 
xCTi)  =  xi 


~  An/.  I  — 


y  d/3 11 


(32) 


where 


fi  =  {(Aq  -  k2i)kziLi  cos  kz  \L[  +  (A02  +  k2x )  sin  kz%\Lx 
-  2A0^,i  sin  A0Ii}/(A02  -  kz2i)2 

+  /{(Aq  -  k2 1 )  A0Li  sin  kz%\Lx  +  2A0Az,i(cos  A0£i  -  cos  frz>1I,)}/(A^  -  ic2i  )2.  (33) 

The  change  in  y  in  the  first  cavity  for  an  electron  with  initial  phase  <f>0  is  given  to  first  order  in  £01 
by 


A-y  __  1 P (£i)  I2  ~  Plo  _  Plo 


yo 


2yo 


yo 


ReUo,/?,*'*0}. 


(34) 


Since  we  have  assumed  that  E(z)  ■=  0  in  the  drift  regions,  we  find  from  Eqs.  (34)  and  (30)  that  the 
transverse  momentum  p  at  the  end  of  the  first  drift  region  is  given  by 


p(L\  +  di,  <f>0)  =  Piol 


1  - 


E0lR{e 


,'*o 


Plo 


/(L|+rf|).VQ  |[0Q— JT]|COS0Q|] 


where 


*„  =  ^^|£0i(7,i  +  *|S,)|, 

yo/3ii 


<Aoi  -  <t>o  +  arg(T}  +  £,</,)  +  arg(E0] )  =  $0  +  §i- 
By  averaging  Eq.  (35)  over  the  initial  phase  4> o,  we  get 


<p(£,  +  «/,)>  -  /]pi0 


/,(*„)  +  i—'R'e —  /0(jrM) 

A. 


-»((Z.,+rf,)A0-8,) 


(35) 

(36) 

(37) 

(38) 


In  writing  Eq.  (38)  we  have  used  the  relations 

~  C  e~ll9~xcos9]  dO  -  i/,Gc), 

2 ir  J0 

J_  f  2"  e*rco*  de  _  y  (x) 

2jr  •'O 


The  presence  of  the  term  A^i  cos  (4>0  +  8])  in  the  phase  in  Eq.  (35)  shows  that  the  electrons 
are  bunched  although  initially  they  have  a  uniform  phase  distribution.  A^i  may  be  interpreted  as 
the  "bunching  parameter"  in  the  first  stage  of  the  amplifier.  If  the  bunching  in  the  first  cavity  is 
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neglected  and  k:  is  set  equal  to  n/L,  then  the  bunching  parameters  Xn  given  by  Eq.  (36)  reduces 
to  the  expression  given  by  Symons  and  Jory  [6}  for  the  two-cavity  gyroklystron.  Although  £0i  is 
small,  we  do  not  make  the  small  argument  expansion  for  J0  Ufn)  and  -M^n)  with  a  view  toward 
the  investigation  of  the  effects  of  nonlinear  inertial  bunching  in  long  drift  regions.  p(L\  +  d\,  <f> o) 
given  by  Eq.  (35)  will  serve  as  initial  condition  for  the  calculation  of  p  and  £  in  the  second  cavity. 


In  the  small  signal  regime,  the  phase-averaged  momentum  ( <p> )  changes  by  a  small 
amount  in  each  cavity  and  the  phase  space  bunching  occurs  mainly  in  the  drift  regions.  Therefore 

in  Eq.  (21),  we  may  make  the  approximation  that  <p>  is  constant  and  replace  it  by  its  value  at 

_  j- 1  _  _ 

the  entrance  of  the  cavity,  i.e.,  at  z  =■  z;  -  £  (L,  +  d,).  With  this  approximation  and  neglecting 

I—  l  __ 

the  homogeneous  solution  (Aj  =  0)  in  Eq.  (21),  we  can  solve  for  p(z)  and  £(z)  by  the  method  of 
successive  approximation  in  all  cavities  j  >  2.  After  lengthy  algebra,  the  following  equations  are 
obtained  correct  to  first  order  in  £0i.  The  transverse  momentum  of  the  electron  at  the  end  of  the 
j'th  stage  is 


P  (f/+l>  *  Pj.o|l 


EniR,e' 


and  the  electric  field  profile  £(z)  in  the  yth  cavity  is 


£(z)  -  E0j  e  ~  ,<V'  G(z,z’)dz\ 


The  complex  quantity  £0/  for  /  ^  2  is  defined  by 

£(,/  *=  -iio  <pG,)> 

-  7oAoU(^i-i)  + 


I  :ty  C-  \  I  <  V  \  EorRr  -  8/.) 

+  '(I  ~  8,,2)  £  — -  e 

I- 1  Plo 


where 


6/  2  =  l  if  /  =»  2 
=  0  if  /  ^  2 

x/*1  =  £  £  Xjy  el(jj\ 
j-  i  r-j 

Xjj  - 


ftoPlO  I  ~ 


I E0lR,(Lj  +  dj)\,  i  <  j 


Xu  =  IXjeitj  -  X,e*‘l,  i  <  j 


£n  =  arg(£01)  +  arg(Tx  +  Rtd  (), 

kjj  “  arg(E0j)  +  arg ( 7}  -I-  Rjdj)  +  8,_ lt  j  >  2 
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*<,  -  arg(E0i)  +  arg(Rj),  j  >  2 

8,- *n 


tan  8j  =  II  XH-  sin  !„•/  it  cos  4V,  y 

1-1 1-1  1-1  i'-l 


01  “  00 


0,  =  (0,  +  Vi)  -  Re  **  I  I  X]T  e‘(jr 


j- 1  >"-/ 


Rj  and  7}  in  Eqs.  (39)  to  (47)  are  given  for  y  >  2  by 


tan  (kz  jLj/lHe^  +  1)  -  —  (e'A°ij  -  1) 

,y  J  An 


*»  -  *i 


r  =  _L  ^  AoA/ tan  (^.iV2)  ~  _  g'A°I;  -  l 

7  A0  A02  -  klj  A02 


2e'A°V2  .  ,  .  A0L,  -  A0Z,;  /c2  ,i 

+  i  ~~2 — rTTT  (Ao  +  ^  sin  -  2A0Arw  cos  — y-  tan  — y 

vA0  —  kz  j)  <■ 


From  Eqs.  (39)  and  (48)  the  terms  Xu  may  be  interpreted  as  generalized  "bunching  parameter." 
Xjj  is  the  contribution  from  self- perturbations  in  the  yth  stage,  and  XtJ  is  the  perturbation  transmit¬ 


ted  from  fth  stage  of  the  amplifier. 


The  average  change  in  y  at  the  end  of  the  yth  cavity  (z  =  £  z,  +  L,)  is  shown  to  be 


<  AX>  -  Re  £f0,  Sk f  Im  £/,(*,_,)  £„,*/• 

y  yo  P\  0  ro 


The  efficiency  rj  =  -  <Ay>/(y0-  l)  in  the  yth  cavity  O'  >  2)  can  be  calculated  by  taking  the 

difference  of  <  -^21  >  between  the  (j  —  l)th  and  the  yth  cavity.  Thus 
y 

V j  =  — (  lm  U, tij-OEojRjl  (52) 

In  Eqs.  (51)  and  (52),  ImGc)  denotes  the  imaginary  part  of  x.  The  output  power  at  the  end  of  yth 
cavity  is 

P out  ”  VjPb’  (53) 

where  Pb  -  VbIb  is  the  beam  power.  The  input  signal  power  is  given  by 

Pm  -  -  ,,P„  (54) 

where  171  is  the  efficiency  in  the  first  cavity,  and  Umn  the  stored  energy  in  TE®„  mode  is 

rr  7T€0m02c5  r^2  _|,  m2 )  ,2/..  ^ 

“  2  ^,lr»\Iw|  2 


/3|2£oi  sin/i  2kzL{  sin  2 kzL{ 
TT  1  2^/Z.,  2k*L\ 


-  -vi 
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where  kx  and  kj  are  the  real  and  imaginary  parts  of  k2  i.  Since  <t>0  is  uniformly  distributed  in  the 
range  0  <  d>0  <  2t r,  the  average  of  Ay  in  the  first  cavity  calculated  from  Eq.  (34)  vanishes.  To 

calculate  17 1,  it  is  therefore  necessary  to  expand  in  the  first  cavity  to  second  order  in  £01,  and 

To 

we  obtain  for  j  =  1 


IK, I 


+  I  mi,  , 


where 


S‘  \AZ-kh2)2  lA°fe-,Ll(  2 k?Lx' 

(Ao  +  k*  i2)£i  f  sinh  2klL\  sin  k?Li 


/' (1  -  cos  2A^|Z,|  1  -  cosh  2kjL\ 

2k*Lx  2  kiLi 


The  amplifier  gain  in  dB  is 


g  =  10  log10 


In  the  small  signal  regime,  both  £out  and  P in  are  proportional  to  £02 1 >  and  g  becomes  independent 
of  /*|N.  From  Eqs.  (53)  and  (57),  it  is  clear  that  PIN  —  0  and  g  —  00  when  1 )]Pb  =  w  UmjQx  with 
m  >  0.  This  corresponds  to  the  self-oscillation  of  the  input  cavity. 

If  the  beam  has  an  axial  velocity  spread,  then  E0j  in  equations  for  p ,  X )j,  r),  and  should 

be  replaced  by  (>3n-av//3M)  E0j.  Furthermore,  an  average  over  the  initial  velocity  distribution  should 

also  be  performed  to  calculate  the  average  values  of  <p>  and  <-^->.  The  expression  for  £0,  in 

y 

Eq.  (41)  will  also  involve  an  average  over  the  initial  electron  velocity  distribution.  In  the  expres¬ 
sion  for  stored  energy  Um  in  Eq.  (46),  /3n  should  be  replaced  by  /3M  av.  The  average  over  the  initial 
distribution  function  is  done  numerically  using  a  Gaussian  distribution. 

We  have  considered  the  interaction  of  the  electron  beam  with  the  TE°mn  modes  of  circular 
waveguide.  However,  the  results  for  TE\\  mode  could  be  applied  to  the  TE^  in  the  rectangular 
waveguide  with  the  following  substitution  [10,11]: 


kn  3f]]/rw  *  ir/ Lx 

c\\k\\Jo(kur^  “ *  -*  / 


cukii'Jo(kuro)  ~~ ’  ~L~L~  Jo(nro/Lx) 

J{(kur0)  — »  J{(nro/Lx ) 

/|(/ciir0)//C|ir0  -*  J^nro/Lj/inro/Lj. 

Lx  and  Ly  are  lengths  of  the  cavity  along  the  x-  and  y-axis. 
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The  small  signal  gain  is  calculated  in  the  next  section  from  Eqs.  (52)  to  (58).  The  difference 
in  phase  between  the  output  and  the  input  signals  can  be  obtained  from  the  phase  of  the  complex 
field  amplitude  E0j  (Eq.  41)  assuming  that  E0l  is  real. 

RESULTS  AND  DISCUSSION 


In  this  section  we  calculate  the  small  signal  performance  characteristics  of  a  three-cavity 
gyroklystron  amplifier  [2]  operating  at  the  fundamental  mode  of  rectangular  cavity.  The 

parameters  of  the  cavities  are:  L\  =  L2  =  0.9x0,  Lj-l.lXo,  d\  —  di  -  i.5\o,  r0-*0.136Xo, 
0i  “  Qi  ”  650,  and  Q 3  =  235.  Lx  of  the  cavities  are  chosen  to  make  the  resonant  frequency  (/<>) 
of  the  cavities  identical  and  Ly/ Lx  =  0.9.  Results  for  a  "cold"  beam  (A  VJ  Vz  =  0)  are  shown  in 
Figs.  2  to  6. 


In  Fig.  2  we  show  the  variation  of  the  small  signal  gain  with  the  magnetic  field  for  three 
values  of  beam  current  (Ib  —  1,  3,  and  6  A).  The  beam  voltage  Vb  is  33.5  kV  and  a  =  v^/v,,  —  1.0. 
The  frequency  /  =  1.0016/o.  As  shown  in  Fig.  1,  amplifier  operation  is  possible  in  different 
ranges  of  the  magnetic  field.  For  stable  amplifier  operation,  the  beam  power,  Pb,  should  be  less 
than  Pbj  (the  threshold  beam  power  for  onset  of  oscillations  in  cavities  /  —  1,2,  3).  In  the  case  of 

11 0 

A  -  6A,  gain  occurs  for  the  magnetic  field  lying  in  the  ranges  5.7  <  2n - <  5.95  and 

ft»0 

IitSIq 

5.98  <  - <  6.3.  The  gain  as  a  function  of  the  magnetic  field  shows  resonance  behavior  at 

COq 


fiekls  where  the  Doppler-shifted  cyclotron  frequency  is  equal  to  the  operating  frequency.  The  peak 
at  fi0  s_5-79  corresponds  to  the  condition  GJ  —  flo/y)L/\n  =  n  for  the  first  two  cavities,  and  the 
peak  at  O0  =  5.95  corresponds  to  this  condition  at  the  last  cavity.  Although  the  gain  is  high  in  this 
magnetic  field  range,  a  stable  amplifier  operation  may  not  be  possible  since  the  gain  is  very  sensi- 

2irfto 

tive  to  variations  in  the  magnetic  field.  For  magnetic  field  in  the  range  5.98  <  -  <  6.3,  the 

a>o 

gain  is  insenstitive  to  variations  of  magnetic  field  and  Pb  <  Pbj  (/  -  1,  2,  3).  Therefore,  stable 
amplifier  operation  with  relatively  high  gain  is  possible  in  this  magnetic  field  range.  In  subsequent 
calculations,  the  parameters  are  optimized  for  maximum  gain  in  the  stable  region  of  operation. 


2nQ0/a)0 


Fig.  2  —  Small  signal  gain  vs  magnetic  field  for 
lb  —  I,  3,  and  6  A.  Other  parameters  are 
y„  -  33.5  kV,  a  -  1.0,  /// o-  1.0016, 

Avz/vr0-  0.0. 


Figure  3  shows  the  small  signal  gain  as  a  function  of  frequency  for  three  values  of 
a  =  vL/V|,  =  1.0,  1.5,  and  2.0  with  Vb  =  33  kV  and  Ib  —  6  A.  The  magnetic  field  is  optimized  for 
maximum  gain  (in  the  region  of  stable  operation)  at  each  a.  The  bandwidth  at  uniform  magnetic 
field  is  extremely  small  (—0.2%).  The  gain  increases  and  the  bandwidth  decreases  as  vx / v(|  is 
increased.  The  magnetic  field  needed  for  maximum  gain  also  increases  with  Vj^/v,,  since  vn 


A. 
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decreases.  The  frequency  of  maximum  gain  /  =  1.0016/0  is  insensitive  to  the  change  in  a.  The 
variation  of  gain  with  a  is  shown  in  Fig.  4  for  /  *  I.OOI6/0,  Vb  —  33.  kV,  and  Ib  —  6  A.  The 
magnetic  field  is  optimized  for  maximum  gain  at  each  a.  The  corresponding  magnetic  fields  are 
also  shown  in  the  figure.  Initially,  the  gain  increases  linearly  with  a  and  then  approaches  satura¬ 
tion.  The  phase  difference  (£)  between  the  input  and  the  output  signals  is  also  calculated  as  a  func¬ 
tion  of  the  various  parameters  such  as  beam  voltage,  current,  magnetic  field,  drive  power,  and  a.  f 
is  found  to  be  a  sensitive  function  of  all  these  parameters  except  the  current.  In  applications  of 
gyroklystrons  requiring  precise  phase  control  such  as  the  RF  linear  accelerator,  it  is  necessary  to 
control  these  parameters  very  carefully.  Figures  5  and  6  show  the  variation  of  £  with  Vb  and 
2rrflo/(oo- 


a 


Fig.  5  —  The  variation  of  the  phase  difference 
(£)  between  the  input  and  output  signals  as  a 
function  of  beam  voltage  for  three  values  of 
magnetic  field.  lb  —  4.0  A,  a  -  1.0,  and 
f/f0-  1.0016. 


11 
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Fig.  6  —  (  as  a  function  of  the  magnetic  field  for 
three  values  of  Vb.  Other  parameters  are  the 
same  as  in  Fig.  S. 


Figures  7  and  8  show  the  effect  of  the  axial  beam  velocity.  A  Gaussian  velocity  distribution 
function  is  assumed,  i.e.,  /(ft,  ft)  <x  exp{(ft  -  ft0)2/2(Aft)2)  8 (ft2  +  pft  -  p$).  The  gain  and  the 
corresponding  optimized  magnetic  fields  are  shown  as  a  function  of  Av2/v,0  in  Fig.  7.  The  gain 
decreases  as  Avz/vr0  increaeses,  and  the  magnetic  field  needed  for  optimum  gain  is  also  increased. 
The  dependence  of  the  small  signal  gain  on  the  magnetic  field  is  shown  in  Fig.  8  for  several  values 
of  AVj/vr0.  As  the  beam  velocity  spread  increases,  there  is  not  only  a  decrease  in  gain  but  the 
range  of  the  magnetic  field  for  stable  amplifier  operation  also  decreases  rapidly.  Stable  operation  of 
a  gyroklystron  amplifier  for  velocity  spread  Av2/v20  >  15%  will  be  difficult. 


Fig.  7  —  Maximum  gain  and  the  corresponding  magnetic  field  as  functions  of  Avz/Vj  for 
a  -  1.0  and  1.5.  Vb  -  33.5  kV,  /,-6A,  and  ///„  -  1.0016 


42 


301 — . — . — . — — . — . — 
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6.3  6.4 


2nS0/«o 


Fig.  8  —  Gain  vs  magnetic  field  for  different  values 
Av/v,  -  0,  0.05,  and  0.10  at  a  -  1.0.  Other  parame¬ 
ters  are  the  same  as  in  Fig.  7. 
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CONCLUSIONS 

We  have  derived  a  comprehensive  small  signal  theory  of  the  multicavity  gyroklystron.  An 
analytic  solution  is  obtained  for  the  "cold"  beam  case.  The  use  of  the  Green's  function  approach 
makes  it  possible  to  satisfy  the  boundary  conditions  for  arbitary  kz  (hence  «),  and  the  gain  as  a 
function  of  frequency  can  be  calculated.  The  stagger-tuned  cavity  configuration  can  also  be  investi¬ 
gated.  The  bandwidth  in  a  uniform  magnetic  field  is  found  to  be  very  small.  The  bandwidth  can 
be  increased  by  using  stagger-tuned  cavities,  but  in  a  uniform  magnetic  field  this  leads  to  a  reduc¬ 
tion  in  gain.  Caplan  [7]  has  shown  that  the  gain  of  the  stagger-tuned  cavity  configuration  is  signifi¬ 
cantly  increased  by  proper  taper  of  the  magnetic  field.  We  plan  to  extend  the  present  theory  to 
include  a  nonuniform  field.  The  present  theory  also  represents  the  first  iteration  of  a  large-signal 
theory  based  on  successive  approximations.  The  small  signal  theory  in  the  bunching  cavities  can  be 
combined  with  a  large  signal  theory  in  the  output  cavity  to  obtain  a  complete  description  of  the 
multicavity  gyroklystron  amplifier.  This  theory  should  become  a  useful  tool  for  obtaining  design 
parameters  of  the  multicavity  gyroklystron. 
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